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The only instance when the General Relativistic (GTR) collapse equations have been solved 
(almost) exactly to explicitly find the metric coefficients is the case of a homogeneous spherical dust 
( Oppenheimer and Snyder in 1939 in Phy. Rev. 56, 455). The equation (37) of their paper showed 
the formation of an event horizon for a collapsing homogeneous dust ball of mass M in that the 
circumference radius the outermost surface, rb ^ ro = IGMji? in a proper time ro oc r'^'^^'^ in the 
limit of large Schrarzschild time t ^ oo. But Eq.(37) was approximated from the Eq. (36) whose 
\mL ■ essential character is (t ~ roln^^y) where, at the boundary of the star y = rs/ro = rhC? /2GM . 

Q,^ , And since the argument of a logarithmic function can not be negative, one must have y > 1 or 

' 2GM/rh(? < 1. This shows that, atleast, in this case (i) trapped surfaces are not formed, (ii) if the 

, collapse indeed proceeds upto r = 0, we must have M = 0, and (iii) proper time taken for collapse 

, ^ . r — > oo. Thus, the gravitational mass of OS black holes are unique and equal to zero. 

O One of the oldest and most fundamental problem of physics and astrophysics is that of gravitational collapse, and, 

. specifically, that of the ultimate fate of a sufficiently massive collapsing body. Most of the astrophysical objects that 
we know of, viz. galaxies, stars. White Dwarfs (WD), Neutron Stars (NS), in a broad sense, result from gravitational 
collapse. It is well known that the essential concept of BHs, in a primitive form, was born in the cradle of Newtonian 
I gravitation In Newtonian gravity, the mass of a collapsing gas cloud is constant even if it is radiating, and equal 
. to its baryonic mass Mi = Mf = Mq. Thus as r decreases, natuarally, the value of M/r = Mq/t steadily increases 
J> ' and at a certain stage one would have 2GM/r = (? when the "escape velocity" from the surface of the fluid becomes 
. equal to the speed of light c (henceforth G = c—\). Soon afterwards, a Newtonian BH would be born. 

■ And in the context of classical General Theory of Relativity (GTR), it is believed that the ultimate fate of sufficiently 
I massive bodies is collapse to a Black Hole (BH). In contrast to Newtonian gravity, here, the gravitational mass of a 

radiating fluid constantly decrease and therefore unlike the Newtonian case, one can not predict with real confidence 
On the actual value of M/ when one would have 2M//r — 1. Given an initial gravitational mass Mi, the three quantities. 
Mi, Mf and Mq are not connected amongst themselves by means of any fundamental constants or by any basic 
I physical principles. Thus, in a strict sense, one reqires to solve the Einstein equations for the collapsing fluid for 
Oh' realistic and ever evolving equation of state (EOS) and radiation transport properties. Unfortunately, the reality is 
Q . that even when one does away with the EOS by assuming the fluid to be a dust whose pressure is zero everywhere 
?H ' including the center, there is no unique solution to the problem. Depending upon the initial conditions like density 
CZJ distributions and assumptions, like self-similarity, adopted, one may find either a BH or a "naked singularity" [^j. 

' It is only when the dust is assumed to be homogeneous, a unique solution can be found Homogeneous or 
^ : inhomogeneous, all dust balls have a very special property: they have no internal energy and they can not radiate. 
Consequuenly, the graviational mass of a dust remains constant during the collapse process, which is, essentially, a 
, Newtonian property. And if the dust collapse starts from a state of infinite dilution at r = r^o = oo, the graviational 

■ mass of a dust must be equal to its baryonic mass M = Mq, another Newtonian property. Having made these remarks, 
we would now proceed to self-consistently analyze the pioneering work of Oppenheimer and Snyder (OS) to determine 
the mass of the BH whose formation it suggested. 

Since the OS solutions are the only (asymptotic and near) exact solutions for the GTR collapse, and are believed 
to explicitly show the formation of an "event horizon" (EH), it is extremely important to critically reexamine them. 
The study of any collapse problem becomes more tractable if one uses the comoving coordinates which are free from 
any kind of "coordinate singularities". By definition, for a given fluid element, the comoving coordinate R is fixed. 
The most natural choice for R is the number of baryons within a given mass shell N{R) or any number proportional 
to it. The comoving time r is of course the time recordrd by a clock attached to a fluid element at R = R. Since a 
dust is in perennial free fall, comoving time is synonymous with "proper time" and, therefore, the dust metric is 

ds^ = dr^ - e'^dR^ - e'^idO^ + sin^ edcp^) (1) 

It turns out that, for any spherically symmetric metric, the angular part is the same and we have 

= (2) 

where r is the invariant circumference radius, and this was the Eq. (27) in the OS paper. OS tried to solve the 



1 



collapse equation by using the metric (2); and by skipping several intermediate equations, we shall focus attention on 
the solutions they oobtained. The Eq. (23-27) of their paper essentially leads to a relationship between r and r: 



2 i?3/2_^3/2 

T 



3 (i?/i?6)3/2r^/2 ' 
By transposing, the foregoing equation yields 



R<Rb (3) 



/ o 1/2 



2/3 



R=['-2^) = ^^^^ 

Apart from the the comoving coordinate system, there is another useful coordinate system, the Schwarzschild system 

ds^ = e^dt^ - e^dr"^ - r^{de^ + sin^ edct)"^) (5) 

and, in principle, it should be possible to work out the collapse (or any) problem in this coordinate system. The 
time t appearing in the Schwarzschild cooordinate system is the proper time of a distant inertial observer. Further, 
by definition, the "comoving" coordinates can not be extended beyond the boundary of the fluid. Thus from either 
consideratiion, the external solutions must be expressed in terms of Schwarzschild system. OS obtained a general 
form of the Schwarzschild metric coefficients which involved derivatives, r = dr/dT and t — dt/dr : 

-5.. = e^ = (l-r')-^ (6) 

and 

gtt = = i-2(i - rf (7) 

Again by skipping few intermediate steps, we note that obtained a general relationship between the coordinate time 
t and coordinate radius R ifor the region inside the collapsing body: 



, 2 /n3/2 3/2 3/2\ o 1/2 I i J/^''^ + 1 

t = i^ro[R^' -r^' y-"' ) -2rQy ' + ro In — (8) 

o y ' — L 

where 



y - I [{R/R^r - 1] + ^ (9) 

It is the above Eq. (8) which corresponds to Eq. (36) in the OS paper, and, for very large t, it attains the form 

Now wc trace some of the intermediate; stcips used by OS and not contained in their paper. By using the simple 
relation ln(m/n) = — ln(n/rn), one may rewrite the above equation as 

yl/2 _ 1 

In the limit of large t and y ^ 1, the above equation becomes 

t = roln(^y^^ (12) 

Or, 

y-l = 46"*/''° (13) 
However, OS overlooked the numerical factor of "4" in their exercise. Prom Eqs. (10) and (14), we find that 
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2 [{R/R^f 



, 3] + ^ = 46-*/- 
roR 



The following equation obtained from the foregoing one is used to eliminate r from relevant equations 

^=4e-*/--^[(i?/i?,)^-3] 



And using Eqs. (15) and (16) into (13) we obtain 

1 
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Rb 
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1 - 



3r, 



1/2, 



2/3 ' 



2R, 



3/2 



(14) 
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However, the corresponding equation in the OS paper (Eq.[37]) contained two small errors: 



-ro In 



Rb 

2n^ 



3r, 



1/2, 



2/3 ' 



2Rl 



(17) 



While the first error is a genuine one; which is due to the omission of the factor 4, the second one is a typographical 
error : the power of in the numerator of the last term of this equation should be 3/2 and not 2. From the foregoing 
equation, they concluded that, "for a fixed value of r as t tends toward infinity, r tends to a finite limit, which increases 
with r". The fact that r is finite at r = ro or r = is evident from Eq. (4) too (provided rb ^ oo and ro > 0). This 
was essentially the idea behind the occurrence of an Event Horizon (EH). 
By differtiating Eq. (16) with respect to r, we obtain 



t = 



(roR\ 
4 \rRbj 



1/2 



Similarly differentiating Eq. (15), and using the above result, we obtain 

. _ R^/yj' 
' Rl'\^/^ 

Now using Eqs. (15), (19) and (20) in Eqs. (7-8), we obtain 



e-^ = 1 - {R/Rbf I 4e-*/'-« + 2 " (^/^&) 



and 



„A-2t/ro 



-t/ro 



+ -[3-(i?/i?,)2] 



(18) 



(19) 



(20) 



(21) 



However, in the paper of OS, these two foregoing equations appear (Eqs. [38-39]) in a slightly erroneous form because 
of the omission of the nemerical factor of 4 in Eq. (13): 



= 1 - [R/Rbf |e-*/'-° + ^ [3 - [R/Rbf] } 



and 



_ -A-2t/ro 



-t/ro 



+ 



l[^-[R/R,f]^ 



(22) 



(23) 



Note these equations were obtained by eliminating r and the t ^ 00 limit covers both the r ^ tq limit as well as the 
further r ^ limit. Now, if we recall that the comoving coordinates R and Rb are fixed, and, when there is a total 
collapse to a physical point at r = 0, the metric coefficients must blow up irrespective of the value of R (i.e, the the 
label of a given mass shell). But this does not happen for the solutions obtained by OS! OS correctly pointed out that. 
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it is only the outer boundary (Rb) for which the (one of the) metric coefficients assume the desired form : ^ oo. 
"For R equal to Rb, tends to infinity like e*/''" as t approaches infinity. But for any interior point, the limitting 
values are different, and infact, = finite even when the collapse is supposed to be complete! OS noted "For R less 
than Rb, tends to a finite limit as t tends to infinity" without trying to find why it is so. 

And on the other hand, when = finite, one sees that e" 0. But, for the external boundary, for which e^ oo, 
it can not be predicted whether e'' indeed approaches zero (these limits remain unchanged even when one ignores the 
numerical factor of 4) . They admitted this problem while they wrote "Also for r < vq, v tends to minus infinity" . 
This 1/ —00 limit would correspond to the singularity i? ^ 0. However they did not ponder why for R = Rb, e" 
does not behave in the desired manner. 

Specifically the fact that for R < Rb, e^ fails to blow up at the singularity definitely hints that there is some tacit 
assumption, made in the OS analysis, which is not realized in Nature ( GTR) or there is a basic fault in the formulation 
of the problem. In our attempt for a possible resolution of this physical anomaly with regard to the unphysical aspect 
of OS solutions, we see from Eq. (8), that, f — > oo if either or both of the two following conditions are satisfied: 

ro ^ 0; i ^ oo (24) 

and 

y ^ 1; t^QO (25) 

OS implicitly assumed that rg is finite and then y ^ 1 and then y < 1 as < ^ oo "For A tends to a finite limit for 
r < ro as t approaches infinity, and for = ro tends to infinity. Also for r < ro, y tends to to minus infinity." But, 
while doing so, they completely overlooked the most important feature of Eq. (8) (their Eq. 36), and Eq. 

(10) that in view of the presence of the t ~ In yi/2li\ term, in order that t is definable at all, one must have 

y > 1 (26) 

For an insight into the problem, we first focus attention on the outermost layer where tjb = Vb/ro, so that the above 
condition becomes 



n > ro (27) 

This condition tells that rb can never plunge below ro- Thus a careful analysis of the GTR homogeneous dust problem 
as enunciated by OS themselves actually tell that trapped surfaces can not be formed even though one is free to chase 
the limit r ^ ro. We may further rewrite this condition as 

!^>1; ^<1 (28) 
ro rb 

This means that, if the collapse indeed proceeds upto rb = 0, the final gravitational mass of the configuration would 
be 



M/(r = 0) = (29) 

But then, for a dust or any adiabatically evolving fluid 

Mi = Mf = constant (30) 

Therefore, we must have Mj = too. This means that ro = in this case, and then, it may be promptly verified 
that, irrespective of the value of r, Eqs. (21) and (22) lead to a unique limiting value for the metric coefficients : 



1 - (46"*/'^" + l) ^ (31) 
or, 

^ oo (32) 

and 

_ gA-2t/ro ^ (33) 
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because I/tq approaches oo much faster than A, if rg = 0. Thus, technically, the final solutions of OS are correct, 
except for the fact they did not organically incorporate the crucial y > condition in the collapse equations 
(here we ignore the missing numerical factor of 4) . And all we have done here is to rectify this colossal lacunae to fix 
the value of ro = 0. 

This conclusion that the work of OS demands — Q could have obtained in a much more direct fashion simply 
from the definition of y in Eq. (15). To this effect, we rewrite this equation as 

y^ha'-l) + — (34) 
2 roa 

Note that during the collapse process R/Rb = a remains fixed for a given mass shell. Suppose we are considering 
the collapse of an interior shell with a < 1 and — 1 < 0. Then if > 0, the second term on the left hand side 
of the foregoing equation can ne made arbitrarily small as r — > 0. Therefore y would become negative as the collapse 
progresses if ro > 0. To avoid this, we must have r/rp > 1. 

We may recall here that, long ago, Harrison et al. [Q (pp. 75) mentioned that spherical gravitational collapse should 
come to a decisive end with Mf = M* = 0. And they termed this understanding as a "Theorem" (without offering a 
real proof): 

"THEOREM 23: Provided that matter does not undergo collapse at the microscopic level at any stage of com- 
pression, then, regardless of all features of the equation of state - there exists for each fixed number of baryons a 
gravitationally collapsed configuration, in which the mass-energy M* as sensed externally is zero." 

In a somewhat more realistic way Zeldovich & Novikov (see pp. 297) discussed the possibility of having an 
ultracompact configuration of degenerate fermions obeying an equation of state p — e/3, where e is the proper 
internal energy density, with A/ — > 0. 

And it is also well known that the so-called "naked singularities" could be of zero-gravitational mass . 

We have already seen that the coordinate time required for collapse to the event horizon or beyond is t = oo. And, 
it may be found that the proper time for collapse of the outer boundary of a dust ball to the central singularity is 0| 

/ r \ 1/2 

r = 7:(^) (35) 
\8mJ ^ ' 

where, the dust ball is assumed to be at "rest" at r = Too at r = i = 0. This equation is also obtainable from Eq. (4) 
provided one chooses i?^ in such a way thay 



RT - ^r^- (36) 



Since M = for the OS problem, the proper time for collapse is infinite: r = oo, and not finite. Physically this means 
that, at any finite proper or coordinate time, there is never any OS-black -hole. The collapse process goes on and on 
indefinitely as spacetime becomes infinitely curved near the would be singularity r = 0. If this picture is correct, for 
self-consistency, Nature should not allow the existence of finite mass BHs. And at this juncture, some readers may 
argue that, this is not so. Irrespective of the results of the OS problem, as highlighted by us in this paper, it may be 
argued that GTR allows the existence of arbitray mass Schwarzschild BHs. 

In another related paper, we assumed the existence of a finite mass BH which is described by Kruskal- Szekeres |^ 
metric. If this assumption of the existence of a finite mass BH is indeed allowed by GTR, one would find the world line 
of a material partcle to be timelike, ds^ > 0, everywhere except probably at the central singularity r = 0. However, 
we have found that, even in the Kruskal- Szekeres metric the metric appears to be null ds^ — 0, at R — 2M. If 
we describe the region interior to the EH by Lemaitre metric ||l0| , we have found that, in this case too, the metric 
becomes null at the EH. We have further shown that, in the Kruskal-Szekeres metric, the metric continues to be null 
for R < 2M too |]. And this can be explained only when we realize that the Event Horizon, itself is the end of 
spacetime for the free falling particle. This implies that mass of the BH is actually M = 0, in complete agreement 
with the present exact analysis. This is also in complete agreement with the idea of Einstein that Schwarzschild 
singularity can not be realized in practice (r = 0) jl^ . 
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